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Lattice QCD simulations have made significant progress in the calculation of nucleon electromagnetic form 
factors in the chiral regime in recent years. With simulation results achieving pion masses of order ~ 180 MeV, 
there is an apparent challenge as to how the physical regime is approached. By using contemporary methods 
in chiral effective field theory (xEFT), both the quark-mass and finite-volume dependence of the isovector 
nucleon magnetic moment are carefully examined. The extrapolation to the physical point yields a result that 
is compatible with experiment, albeit with a combined statistical and systematic uncertainty of ±10%. The 
extrapolation shows a strong finite-volume dependence; lattice sizes of L > 5 fm must be used to simulate 
results within 2% of the infinite-volume result for the magnetic moment at the physical pion mass. 
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I. INTRODUCTION 

The distribution of the electric and magnetic charge cur- 
rents of the nucleon are characterized by the elastic electro- 
magnetic form factors — for recent reviews on experimen- 
tal progress, see Refs. fl-Ht]. The description of the electro- 
magnetic form factors in terms of the elementary degrees of 
freedom of QCD has seen significant progress through recent 
advances in lattice QCD simulations flq- UOIl . 

Lattice QCD simulations of the electromagnetic form fac- 
tors of the nucleon 111 111 are now probing into the chiral regime, 
where the QCDSF Collaboration have recently reported re- 
sults at pion masses as low as ~ 180 MeV iflOll . The re- 
sults of this work have presented a challenge in the pion-mass 
extrapolation to the physical point. The results presented in 
Ref. lUOll have been used to investigate the applicability of a 
range of chiral effective field theory (xEFT) methods, includ- 
ing "heavy-baryon" [12], "small scale expansion" lfl3ll . and 
"covariant baryon" approaches It has been demonstrated 
that there is a difficulty in achieving a consistent quantitative 
description of the pion-mass dependence of key observables 
(such as magnetic moments and charge radii) between the 
physical point and the lightest simulation results [10]. This 
issue has persisted across a number of simulations Il5lll6ll . In 
the present manuscript, a chiral extrapolation of the QCDSF 
results is developed for the isovector nucleon magnetic mo- 
ment based on finite-range reg ularized (FRR) vEFT JTtL fl8h . 

In the application of FRR to the extrapolation of the mag- 
netic moment, new developments are utilized in order to en- 
sure the robustness of the extrapolation procedure. These in- 
clude the identification of the preferred finite regularization 
scale directly from the lattice results lfl9ll and a determination 
of an upper bound of the pion mass that can be reliably in- 
corporated in the extrapolation [20]. A feature of the analysis 
is that the combination of both finite-volume corrections and 
the onset of rapid nonanalytic behaviour in the chiral regime 
leads to an extrapolation that is compatible with the experi- 
mental value. The analysis provides the predicted quark-mass 
dependence for a range of fixed-volume lattices, which act to 
emphasize the importance of achieving large volumes in or- 



der to reveal the strong nonanalytic behaviour directly on the 
lattice. 

The lattice QCD results for the magnetic moment of the 
isovector nucleon, fi^, from the QCDSF Collaboration are 
displayed in Fig. Q] The lattice calculation used Nf = 2 and 
the 0(a)-improved Wilson quark action II X Oil . The isovector 
combination (p — n) is considered to avoid calculating the 
disconnected loops that occur in full QCD. To ensure that the 
lattice results give a reasonable approximation to the infinite- 
volume limit, the following restrictions are applied: L > 1.5 
fm and m^L > 3. There are nine lattice points that satisfy 
these criteria from the original set of results. The lattice sizes 
considered vary from 1.7 fm to 2.9 fm. A simple linear fit 
is included in this plot, which does not take into account the 
chiral loop integrals, nor the finite-volume corrections to the 
data. Neglecting these effects, it is not surprising that the lin- 
ear trend does not reach the experimental value of the mag- 
netic moment at the physical pion mass. The use of extended 
XEFT methods in performing the extrapolation to the physical 
point will now be explored. 
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FIG. 1. (color online). Lattice QCD data for fi v N from QCDSF |l(J], with 
experimental value as marked [21, 22]. The lattice results satisfy L > 1.5 fm, 
m w L > 3. A simple linear fit is also included, which misses the experimental 
value. 
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II. CHIRAL EFFECTIVE FIELD THEORY 

The elastic matrix element for the baryon-photon interac- 
tion can be parametrized by the Dirac and Pauli form factors, 
F\ and F 2 , respectively, written as: 

{B(j>')\J,\B{p)) = 

u s \p') {7, F,{Q 2 ) + F 2 {Q 2 )\ u'(p) . (1) 



This expansion contains an analytic polynomial in m 2 plus 
the leading-order chiral loop integrals (Tn,a), from which 
nonanalytic behaviour arises. The coefficients af are the 
(scale-dependent) 'residual series' coefficients. Upon renor- 
malization of the divergent loop integrals, these terms corre- 
spond to low-energy coefficients of %EFT [13 ill . In this in- 
stance, only two free parameters are provided in the residual 
series. The leading-order diagrams included in this investiga- 
tion are simply the 1-meson loops, as shown in Figs. |2]and|3] 



Q 2 is a positive momentum transfer Q 2 = —(p'—p) 2 . The 
Sachs electromagnetic form factors Ge,m are the linear com- 
binations of F% and F2 defined by: 



G E (Q 2 ) = FxiQ 2 ) 
G M (Q 2 ) = ^i(Q 2 ) 



Q 2 

4m 2 B 

F 2 {Q 2 ). 



F 2 (Q 2 ), 



(2) 
(3) 



The Sachs magnetic form factor of the nucleon at zero- 
momentum transfer, Gm{Q 2 — 0), defines the magnetic mo- 
ment as two separate terms: the Dirac moment (unity), plus an 
anomalous contribution associated with the internal structure 
of the hadron: 



= 1 + K n . 



(4) 
(5) 



For the leading-order contributions to the magnetic mo- 
ment, the standard first-order interaction Lagrangian from 
heavy-baryon chiral perturbation theory (%PT) is used 1I23I - 
IH: 

4pt = 2DTt [B v S^{A^, B v } ] + 2FTr [B v Si^[A^ B v ] ] 
+ C(f^A^B v + B v A^), (6) 

(7) 



where the pseudo-Goldstone fields are encoded as the adjoint 
representation of SU(3) l ® SU(3)a, forming an axial vector 
combination A,, : 



(8) 
(9) 



£ = exp 



By the convention presented here, f n = 92.4 MeV. The values 
for the couplings in the interaction Lagrangian are obtained 
from the SU(6) flavour-symmetry relations ll25ll29ll and from 
phenomenology: D = 0.76, F = |D and C = -2D. 

From the full Lagrangian, the chiral behaviour of the mag- 
netic moment can be written in terms of an ordered ex- 
pansion in pion mass squared, through use of the Gell- 
Mann— Oakes— Renner Relation, m q oc m 2 , O30ll : 



a£ + 4 ml + T N {ml ; A) + T^m 2 ; A) + 0{m%) . 

(10) 



A. Finite-range regularization 

Finite-range regularized effective field theory handles di- 
vergences in the ultraviolet region of the loop integrals by in- 
troducing a regulator function as part of the coupling to each 
vertex of the diagram. The regulator function u(k ; A) intro- 
duces a cutoff scale A, and should satisfy u|fe = o = 1 and 
w|fc->oo = 0- Th e exact functional form chosen for the regu- 
lator is independent of the result of calculation, so long as the 
lattice QCD results are constrained within the power-counting 
regime (PCR). A smooth attenuating dipole form is chosen for 
this investigation: 



u(k;A) 



1 



A2 



(11) 



Detailed analyses exist for a range of alternative forms lfl9[ 
l32ll . Though it has been suggested that a sharp cutoff FRR 
scheme should be chosen to ensure the preservation of the chi- 
ral Ward Identities l33ll . it is possible to maintain chiral sym- 
metry by including the necessary vertex corrections at higher- 
order in order to suppress any extra scheme-dependent non- 
analytic terms induced by regulators such as the dipole lf34ll . 
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FIG. 2. The pion/kaon loop contributions to the magnetic moment of an octet 
baryon B, allowing a transition to a baryon B' , with a photon attachment, 
which provides the leading nonanalytic contribution. All charge conserving 
transitions are implicit. 
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FIG. 3. The pion/kaon loop contribution to the magnetic moment of an 
octet baryon B, allowing transitions to nearby and strongly coupled decuplet 
baryons T. 
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Here, chiral Ward Identities are maintained to the order of the 
calculation. 

FRR, in conjunction with xEFT-inspired techniques, pro- 
vides a robust method for achieving an extrapolation to physi- 
cal quark masses, and identifying an intrinsic scale embedded 
within lattice QCD results. It has been demonstrated previ- 
ously that such an intrinsic scale may be extracted from the 
results of lattice QCD calculations for the mass of the nucleon 
lfl9ll . This property is a consequence of the size of the PCR, 
defined where the expansion formulae of %PT formally hold 
(to finite chiral order). This extended effective field theory 
proceeds by analyzing the behaviour of the renormalization 
of one or more low-energy coefficients of the chiral expansion 
as a function of the regularization scale. Ideally, with lattice 
QCD results constrained entirely within the PCR, the renor- 
malized coefficients are independent of regularization scale. 
However, in practice, a scale-dependence is observed; par- 
ticularly for lattice result sets including points corresponding 
to quark masses beyond the PCR. By truncating the lattice 
QCD results at different quark-mass points (corresponding to 
a value of to 2 max ), an optimal regularization scale can be 
identified. This optimal scale is the value at which the low- 
energy coefficients are least sensitive to the truncation of the 
lattice results fijj||. 



B. Loop integrals and definitions 

The value of the magnetic moment is renormalized by con- 
tributions from loop integrals, obtained from the effective field 
theory. Here, the focus is on the pion contributions, noting 
that it is straightforward to include the kaon contribution, us- 
ing m 2 K = m 2 K phys + § (m 2 . - «4, phys ), as is done in Sec. HI 
reporting the results. The two leading-order loop integrals are 
the 1-meson loops, as shown in Figs. [2]and[3] Each loop inte- 
gral can be expanded as a polynomial series, which is analytic 
in quark mass, plus a nonanalytic term: 



T N (ml;A) 
r A {ml;A) 



0(m\ 



(12) 



bo + Xn + b% n 

+ b% ml + X A mi log m^/fj, + 0(m%) 



(13) 



Here, \x is an implicit mass scale, chosen here to be 1 GeV. The 
coefficients of each polynomial, bf' A are entirely dependent 
on the choice of finite-range regularization scheme, and so 
they are regulator-dependent quantities. The renormalization 
program of FRR combines the scale-dependent b{ coefficients 
from the chiral loops with the scale-dependent a, coefficients 
from the residual series in Eq. ( [Tol l at each chiral order i, such 
that the result is a scale-independent coefficient 



Co 

C2 



UN 



UN 



etc. 



(14) 
(15) 



This means the underlying a* coefficients undergo a renor- 
malization from the chiral loop integrals. The renormalized 
coefficients Cj are an important part of the extrapolation tech- 
nique. A stable and robust determination of these parameters 
forms the heart of determining an optimal scale A scale . 



The loop integrals can be expressed in a convenient form by 
taking the nonrelativistic limit and performing the pole inte- 
gration for fco- Renormalization is achieved by subtracting the 
relevant £>q term from the integral, effectively absorbing it into 
the correspon ding renormalized coefficient cq. The integrals 
take the form [28, 351: 



T N (ml;A) 
7a (m 2 ; A) 



~XN 

3tt 2 

-XA 

3tt 2 



d 3 fc fc 2 U 2 (fc;A) 



(16) 



d*k 



(k 2 + m 2 ) 

k 2 {2w{k) + A)u 2 (k-A) A 

2w 3 (fc) [w(fc) + A] 2 ' 
(17) 



The chiral coefficients xn and xa are determined from inter- 
actions in the chiral Lagrangian of Eq. ([6}: 



p 

Xn = 



P 



m N 



8^ 
mjv 2C 2 



(D + F) 2 = - Xn , 



-Xa- 



(18) 



(19) 



The chiral expansion of the magnetic moment in Eq. [10] can 
now be written out in a form renormalized to order 0(1): 

= c + 4 ml + t N (ml ; A) + 7a (m 2 . ; A) + 0(m%) . 

(20) 

Since lattice simulations are necessarily carried out on a 
discrete spacetime, any extrapolations performed should take 
into account finite-volume effects. ;\EFT is ideally suited for 
characterizing the leading infrared effects associated with the 
finite volume. In order to achieve this, each of the three- 
dimensional integrals can be transformed to its form on the 
lattice using a finite-sum of discretized momenta, see Allton 
etal. 0611 for instance: 

(2nf 



E 



L x L y L z , 

"'X ; 1/ ' "'2 



(21) 



On the finite-volume lattice, each momentum component is 
quantized in units of 2ir/L, that is fe, = 2irrii/L for integers 
m. Finite-volume corrections S FVC can be written simply as 
the difference between the finite sum and the corresponding 
integral. It is known that the finite-volume corrections satu- 
rate to a fixed result for large values of regularization scale 
Clllll. The value A' = 2.0 GeV is chosen to evaluate all 
finite-volume corrections independently from the integral cut- 
off scale A in Eqs. ( TTol ) and ( fTTI i. This method is equivalent to 
the more algebraic approach outlined in Ref. [38]. The finite- 
volume version of Eq. ( f20b can thus be expressed: 

Hi = co + 4 ml + (f N (ml ■ A) + A')) 
+ (f A (m 2 ; A) + 5™ G {ml; A')) + 0{m\) . (22) 



III. RESULTS 
A. Renormalization flow analysis 

In order to obtain the most robust extrapolation, an optimal 
regularization scale is sought. The robustness of the extrap- 
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olation is characterized by its stability against truncation of 
the lattice data set. That is, a similar extrapolation should be 
achieved regardless of the number of data points used in the 
fit. The optimal regularization scale A scal ° may be obtained by 
calculating the low-energy coefficients (e.g. Co) from Eq. (l22l 
for a range of regulator values A. Since the lattice simulation 
results extend outside PCR, the renormalized value of the co- 
efficients will be scale dependent. However, the analysis in 
Ref. 11911 demonstrates that using different amounts of lattice 
data yields a different scale dependence. If the lattice simula- 
tion results lie close to the PCR, the scale dependence is nat- 
urally less than if a more extensive set of lattice data is used. 
The optimal regularization scale is the value of A at which the 
same value of cq (or any low-energy coefficient) is obtained 
regardless of the amount of lattice data used. This is the scale 
where the values of the low-energy coefficients correspond to 
the values obtained from lattice results within the PCR I19I1 . 

Consider the behavior of cq from Eq. ( l22l i as a function 
of the regularization scale A. Using different upper values 
of max , a set of renormalization flow curves may be con- 
structed. The renormalization flow curves including up to all 
nine lattice results, and using a dipole regulator, are plotted on 
the same set of axes in Fig. [4] As more data are included in 
the fit, a greater degree of regulator dependence is observed. 
Note that there is a reasonably well-defined A value at which 
the renormalization of cq is least sensitive to the truncation of 
the data. This indicates that there exists an optimal regulariza- 
tion scale embedded in the lattice QCD results. 

The optimal regularization scale for a dipole can be ex- 
tracted from Fig. [4] using a x\ f analysis. Such an analysis 
will also provide a measure of the systematic uncertainty in 
the optimal regularization scale. By plotting Xdof a g amst A, 
where dof equals the number of curves n minus one, a mea- 
sure of the spread of the renormalization flow curves can be 
calculated, and the intersection point obtained. The x\ Q f * s 
constructed at each value of A for cq (with uncertainty Scq): 



Xdof 



i 



is 

i=l 



(4(A)-cp(A)r 
(<H(A)) 2 



co (A) 



£r = i4(A)/(*4(A)r 
e;u v(5c (a)) 2 



(23) 
(24) 



The indices i and j correspond to data sets with different val- 
ues of max . The Xdof Pl Qt corresponding to Fig. |4]is 
shown in Fig. [5] The upper and lower bounds on A obey the 
condition x 2 dof < X d of,min + 1 /( do f)- Thus the optimal reg- 
ularization scale for a dipole form is: A scalc = 0.82+q']J GeV. 
This value is smaller than the optimal regularization scale ob- 
tained for the nucleon mass using a dipole form |19] based 
on lattice QCD results from JLQCD & PACS-CSffl and 
CP-PACS [4111 . Nevertheless, the value of the optimal scale 
is the same order of magnitude as those calculated from the 
nucleon mass analyses, and previous studies of the magnetic 
moment 14211 . In addition, it has been assumed that all avail- 
able lattice results should be used in constructing the fit. This 
is not necessarily the case. The method described in Ref. ll20ll 
outlines a procedure for determining the optimal fit window 



of pion masses, as discussed in the following section. 



B. Chiral extrapolations 

Using the optimal regularization scale, extrapolations can 
be made on various lattice volumes. In order to determine the 
most suitable number of data points to be used for fitting the 
lattice results, one may perform several extrapolations, using 
a varying number of points each time, and compare the size 
of the statistical and systematic error estimates. The statistical 
uncertainty comprises contributions from the fit coefficients, 
and the optimal regularization scale A scalc , which is assumed 
to be independent of the other fit coefficients, and its contri- 
bution is added in quadrature. The axial coupling and the 
pion decay constant f v are assumed to be sufficiently well de- 
termined experimentally. The systematic uncertainty in the 
extrapolation is estimated by comparing the results from dif- 
ferent regulator functional forms. The triple-dipole regulator 
introduced in Refs. Ill 9U20I1 is considered, which interpolates 
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FIG. 4. (color online). The renormalization flow of co for fj?L obtained 
using a dipole regulator on QCDSF lattice QCD results. For each curve, two 
arbitrary values of A are chosen to indicate the general size of the error bars. 
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obtained using a dipole regulator on QCDSF lattice QCD results, up to and 
including ra\ max = 0.81 GeV 2 . 
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FIG. 6. (color online). Behaviour of the extrapolation of fj,^ to the physical 
point vs m 2 max . In each case, the value of A acale is used, as obtained from 
the corresponding x\ a f analysis. The error bars include the statistical and 
systematic uncertainties added in quadrature. 

between the dipole and the sharp cutoff regulators. The resul- 
tant chiral extrapolation using a dipole regulator is compared 
to that using a triple-dipole regulator in order to estimate the 
systematic uncertainty. 

The quadrature sum of the statistical and systematic uncer- 
tainties in the extrapolation of n v N to the physical point, for 
different values of m 2 max , is shown in Fig. [6] Ideally, one 
should find a best value of the upper limit to 2 max , as indi- 
cated by the best compromise between statistical and system- 
atic effects. Fig. [6] indicates that the smallest error bar oc- 
curs when all nine lattice points are included. This set cor- 
responds to a maximum pion-mass value: m 2 max = 0.81 
GeV 2 . However, it is helpful to know the relative contribu- 
tions from statistical and systematic sources. Table H] summa- 
rizes the breakdown of each error bar into its sources. Fig. 
[7] shows the magnitude of the statistical and systematic error 
bars for different values of m 2 max . Clearly, the general trend 
of the statistical error bar decreases as more lattice results are 
considered, and likewise the general trend of the systematic 
error bar increases. At m 2 _„ v = 0.44 GeV 2 , the statistical 
and systematic error bars are closest in magnitude, which indi- 
cates the proximity of a 'sweet spot' (denoted m 2 max = m 2 ), 
at which the best trade-off between statistical and systematic 
uncertainty is achieved. 

The renormalization flow curves corresponding at most to 
a value of m 2 max = 0.44 GeV 2 are shown in Fig. [8] and 

the corresponding Xdof P^ oi * s snown m Fig- Clearly, the 
statistical contribution to the uncertainty in the extrapolation 
is larger than in the case where all lattice results are used. 
This is reflected in the larger uncertainty in the identification 
of the optimal regularization scale (at optimal to 2 max = m 2 : 
Ascf le _ 0.87^036- This value is consistent with the opti- 
mal regularization scale obtained for the nucleon mass, using 
a dipole form 01911 . This provides evidence for the successful 
extraction of the intrinsic scale in the nucleon-pion interac- 
tion. Since this value of A s _ c t lc is obtained from the best com- 

m 

promise between statistical and systematic uncertainty, it will 
be used in the following chiral extrapolation of the isovec- 



TABLE I. Results for the isovector nucleon magnetic moment for different 
values of max , extrapolated to the physical point, corresponding to Fig. 
[6] The uncertainty in fj,^ (m^ phys ) is provided in the following order: the 

statistical uncertainty, the optimal regularization scale A scalc , and the sys- 
tematic uncertainty due to the regulator functional form, respectively. 

"4,max(GeV 2 ) MivO^.phys) OmO 

0.185 4.35(13) (57) (17) 

0.230 3.79(10)(60)(33) 

0.240 3.79(9)(61)(32) 

0.436 4.29(7) (34) (40) 

0.449 4.29(7) (29) (40) 

0.810 4.26(5)(16)(39) 
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FIG. 7. (color online). Magnitude of the statistical and systematic error 
bar in the extrapolation of to the physical point vs m 2 max . In each 
case, the value of A scalc is used, as obtained from the corresponding xjj / 
analysis. At a maximum pion mass of m 2 = m 2 max = 0.44 GeV 2 , the 
best compromise between statistical and systematic uncertainty is achieved. 

tor nucleon magnetic moment. Only the lightest seven lattice 
points (corresponding to to 2 max = 0.44 GeV 2 ) will be used 
in the fit. 

Consider the behaviour of the magnetic moment as a func- 
tion of the quark mass. Extrapolation curves corresponding to 
infinite volume, and a variety of finite volumes are shown in 
Fig. [10] For each curve, only the values for which m^L > 3 
are plotted. These finite volumes include typical values at 
which current lattice QCD results are produced. For example, 
a full QCD simulation at physical quark masses on a (4 fm) 3 
volume will significantly underestimate the infinite-volume 
result. These curves indicate that a box length of L > 5 fm is 
required to achieve an extrapolation within 2% of the infinite- 
volume result. 

The finite-volume expansion of Eq. ( 1221 is constrained by 
the lattice simulation results from several different volumes 
in the range 1.7 — 2.9 fm, as shown in Fig. [TJJ The infinite- 
volume extrapolation is fit to the lattice simulation results only 
after the results have been corrected to infinite volume. These 
points are shown in Fig. [T_2] The extrapolation to the physical 
point also includes an inner error bar representing only the sta- 
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FIG. 8. (color online). The renormalization flow of co for fi v N obtained 
using a dipole regulator on QCDSF lattice QCD results. Only the lightest 
seven lattice results are used. For each curve, two arbitrary values of A are 
chosen to indicate the general size of the error bars. 
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lattice sizes: 1.7 — 2.9 fm. The experimental value is marked [21, 22|]. In all 
finite-volume extrapolations, the provisional constraint m^L > 3 is used. 
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tistical uncertainty, and an outer error bar, which also includes 
the systematic uncertainty due to the regulator in quadrature. 

In all extrapolations, the strange quark loops have been un- 
quenched, and the effects of kaons loops that would occur 
in an SU(3) lattice calculation are estimated. The result is 
a change of only w 0.7% larger at the physical point when 
kaons loops are included. 

The finite-volume extrapolations of Fig. Q7j]are generally 
useful for estimating the result of a lattice QCD calculation 
at certain box sizes. This can also provide a benchmark for 
estimating the outcome of a lattice QCD simulation at larger 
and untested box sizes. 
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FIG. 11. (color online). Extrapolations of fi^ at different finite volumes 
and infinite volume. The finite-volume lattice QCD results from Ref. floll 
are plotted for comparison, with box sizes in the range 1.7 — 2.9 fm. In all 
finite-volume extrapolations, the provisional constraint m, L > 3 is used. 
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FIG. 12. (color online). Extrapolations of [i v N at different finite volumes and 
infinite volume. The lattice QCD results displayed have been corrected to infi- 
nite volume. Only the lightest seven points are used in the fit, corresponding 
0.44 GeV 2 . 
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IV. CONCLUSION 

The technique for obtaining an optimal regularization scale 
from lattice QCD results has been investigated in the context 
of the magnetic moment of the isovector nucleon, using recent 
precision lattice QCD results from QCDSF. An optimal reg- 
ularization scale was identified by analyzing the renormaliza- 
tion flow of the low-energy coefficient cq with respect to the 
scale A, whilst extending beyond the power-counting regime. 
An optimal value of max was also obtained, where the sta- 
tistical and systematic error estimates of a chiral extrapolation 
are comparable in magnitude. This value rh 2 provides a guide 
to the range of pion masses in which finite-range regulariza- 
tion techniques are not the dominant source of uncertainty in 
a chiral extrapolation. 

A regularization scale A^§ lc was determined where the 
renormalization of Co is least sensitive to the truncation of the 
lattice QCD results. The value of the optimal regularization 
scale was consistent with results from the nucleon mass anal- 



ysis. Thus an intrinsic scale has been uncovered, which char- 
acterizes the energy scale of the nucleon-pion interaction. The 
result therefore further demonstrates the success of the proce- 
dure for using lattice QCD results to extrapolate an observable 
to the low-energy region of QCD. 

Using the value of the intrinsic scale, the extrapolation of 
the magnetic moment to the physical pion mass and infinite- 
volume lattice box size is consistent with the experimental 
value. More importantly, the finite-volume extrapolations pro- 
vide a benchmark for estimating the outcome of a lattice QCD 
simulation at realistic or optimistic lattice sizes. This serves 
to emphasize the importance of achieving large volumes in 
realizing the correct nonanalytic behaviour found in nature. 
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